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Abstract. To addres s the important issue of how kinetic energy of collimated blast waves is converted into radiation, Pohl 
& Schlickeiser (2000) have recently investigated the relativistic two-stream instability of electromagnetic turbulence. They 
have shown that swept-up matter is quickly isotropized in the blast wave, which provides relativistic particles and, as a result, 
radiation. Here we present new calculations for the electrostatic instability in such systems. It is shown that the electrostatic 
instability is faster than the electromagnetic instability for highly relativistic beams. However, even after relaxation of the beam 
via the faster electrostatic turbulence, the beam is still unstable with respect to the electromagnetic waves, thus providing the 
isotropization required for efficient production of radiation. While the emission spectra in the model of Pohl and Schlickeiser 
have to be modified, the basic characteristics persist. 
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1. Introduction 

Published models for the 7-ray emission of blazars are usually 
based on interactions of highly relativistic electrons of unspec- 
ified origin propagating in the relativistic jet of active galactic 
nuclei (Hartman et al. 1997, Bicknell et al. 2001, Catanese and 



structures which relativistically move along a common trajec- 
tory and can be followed over years. On larger scales optical 
data still show individual knots, but also emission in-between 



Weekes 1999). The usual shock or stochastic electron accelera- 



tion processes would have to be very fast to compete efficiently 
with the strong radiative losses at high electron energies, which 
operate on timescales of the order of ten seconds for the rapidly 
variable 7-ray blazars, which is similar to the duration of a sin- 
gle shock crossing cycle for TeV electrons at a moderately rela- 
tivistic shock assuming Bohm diffusion (Gallant & Achterberg 



1998). Energetic electrons may also result as secondaries from 



photomeson production of highly relativistic hadrons, but the 
observed short variability timescale places extreme conditions 
on the magnetic field strength, for the hadron gyroradius has 
to be much smaller than the system size, and on the Doppler 
factor, for the intrinsic timescale for switching off the cascade 
is linked to the observed soft photon flux. 

Here we consider the energisation of relativistic particles 
in jets by interactions with the surrounding medium. The rapid 
variability displayed in the 7-ray light curves of blazars re- 
quires the emission regions be less than ^0.01 pc in size. 
Also, VLBI observations of blazars indicate that the jets are not 
continuously filled emission regions, but consist of individual 
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(e.g. Boksenberg et al. 1992 ). Detaile d spec tral studies (e.g. 
Meisenheimer, Roser and Schlotelburg 1996 ) indicate that the 
spectrum of the synchrotron radiating electrons is fairly inde- 
pendent of position, suggesting that many acceleration sites 
with uniform characteristics are operational on scales of ~ 100 
pc. 

Accounting for these findings we model a jet as a chan- 
neled outflow with relativistic bulk velocity V, consisting of 
isolated plasma clouds or blast waves which contain cold elec- 
trons and protons of density ri}, (see Fig.[j]). These clouds may 
correspond to the individual components observed in VLBI im- 
ages of blazars, whereas many of them and secondary electrons 
of escaped nucleons would contribute to the optical and X-ray 
jets observed on larger scales in nearby AGN such as M87. 
For convenience we assume that the outflow is directed paral- 
lel to the uniform background magnetic field. We also neglect 
a possible expansion of the jet cloud in this paper. Expansion 
is expected to occur on much longer timescales than the effects 
discussed here. 

Each cloud of protons and electrons propagates into the 
surrounding interstellar medium that consists of cold protons 
and electrons of density n* (quantities indexed with an aster- 
isk are measured in the laboratory (galaxy) frame). Viewed 
from the coordinate system comoving with the outflow, the 
interstellar protons and electrons represent a proton-electron 
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Fig. 1. Sketch of the basic geometry. The thickness of the chan- 
neled blast wave d, measured in its rest frame, is much smaller 
than its halfdiameter. The blast wave moves with a bulk Lorentz 
factor r through ambient matter of density rii. 

beam propagating with relativistic speed —V antiparallel to 
the uniform magnetic field direction. This situation is unsta- 
ble and waves will be excited which backreact on the incoming 
beam. In our earlier analysis (Pohl & Schlickeiser 2000 1 the 
stability of this beam was examined under the assumption that 
the background magnetic field is uniform and directed paral- 
lel to the direction of motion. It was shown that the beam very 
quickly excites low-frequency electromagnetic waves, which 
quasi-linearly isotropize the incoming interstellar electrons and 
protons in the blast wave plasma, thus providing relativistic 
particles. These processes are essentially those thought to oper- 
ate in the first half-cycle of relativistic shock acceleration, but 
rather than assuming the existence of a hydrodynamical shock 
and a particular intensity of MHD turbulence, kinetic theory 
was used to self-consistently calculate the build-up of turbu- 
lence and its backreaction on the incoming interstellar parti- 
cles. 

The protons carry the bulk of the power, because veloci- 
ties are conserved in the isotropization process and hence the 
protons have a factor of ~ 2000 more kinetic energy than the 
electrons. The protons produce high energy emission via the 
production of neutral pions in inelastic pp collisions. In par- 
allel charged pions would provide secondary electrons which 
emit synchrotron radiation as well as bremsstrahlung and in- 
verse Compton scattering. As discussed in Pohl & Schlickeiser 
( 2000 ), the high energy emission produced by these parti- 
cles has characteristics typical of BL Lacertae objects. The ob- 
served secular variability of the 7-ray emission of AGN would 
be related to the existence (or non-existence) of a relativistic 
blast wave in the sources, and thus to the availability of free 



energy in the system. The observed fast variability, on the other 
hand, would be caused by density inhomogeneities in the inter- 
stellar medium through which the blast wave propagates. 

In parallel to the 7-rays , neutrinos are emitted whose spec- 
trum and flux would be closely correlated with those of the 7- 
rays, which permits one to use the 7-ray light curves of blazars 
to very efficiently search for neutrino emission as a diagnostic 
for an hadronic origin of the high energy radiation (Schuster et 
al. pOOll). 



In the initial stage virtually no radio emission would be 
produced, for the plasma frequency and the synchrotron self- 
absorption frequency would be too high. After the initial 7-ray 
emission phase, when the blast wave has decelerated and ex- 
panded, a mm-radio flare would build up, which behaviour is 
preferentially observed in radi o to 7 -ray correlation studies of 
EGRET sources (Miicke et al. 



1996) 



Here we expand on the previous treatment by calculating 
the two-stream instability for longitudinal, electrostatic waves. 
In contrast to electromagnetic waves, which scatter the parti- 
cles in pitch angle but preserve their kinetic energy until the 
distribution is isotropized, the electrostatic waves change the 
particles' energy until a plateau distribution is established. 

2. Longitudinal instabilities of a proton-electron 
beam 

2. 1. The dispersion relation 

For charge e, mass m particles under the action of an electrical 
field E = eii-Eii exp(tk(x — at)), in the direction parallel to 
an ambient magnetic field the perturbations, 5f, to the original 
distribution function satisfy 



dSf dSf 



dt 



dx 



eE 



dp 







(1) 



which, to first order in <5///o, and with Sf also taken to vary as 

5f(p) exp(tk(x — at)) gives 



ifc(«i| -a)*/ + e£7|i|^ = 



i.e 



dp\\ 
ieEu 



df 



k(v\\ ~ a) dp\\ 
The electrostatic balance equation requires 

div.E = 47r^e J Sfd 3 p 

where the sum is over charge species. Thus 



fc 2 £|| = 47r£|| J2< 



dfo 

9P " d 3 p 



u„ -a 

Then the dispersion relation is 

9/o 



(2) 



(3) 



(4) 



(5) 



■d 3 p 



(6) 
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Let a = a R + iaj. Then split the RHS of eq. |] formally into the 
real part, J, and the imaginary part, /, 



k 2 = J(a R , ai) + iI(a R , a\) 
Because k 2 is meromorphic, then 

dI(a K ,ai) 



I(a R , a{) = ~ I(a R , 0) + ai- 

The Cauchy condition yields 

dI(a R ,ai) _ dJ(a R ,ai) 
dai da R 

Thus the lowest order approximation yields 



ai=0 



«I 



-/(a R ,0) 



( dJ(a R ,0) 



Then 

J(a R ,0) = 4vr^e 2 P 



Op i 



d 3 p 



(7) 



(8) 



(9) 



(10) 



(11) 



v\\ - a R 

where P denotes the Cauchy principal value, while 

I(a R , 0) = 4^ 2 J> 2 S ( v \\ - a «) d3 P (l V 



Using the notation of Pohl & Schlickeiser ( |2000t ) the initial dis- 
tribution function in the blast wave frame is (after correction of 
a typo in their Eq.3) 



f(p,t = 0) 



riiS(px)S(p\\+P) n h S(p±)5(pi\) 



2np± 

In this case Eq. ^[reduces to 



27rpj 



k z = Q 2 



S 



(13) 



(14) 



: (a + V) 2 _ 
where Q 2 = u 2 p + uj 2 c and 5 = m/ (n b T 3 ). 

2.2. Properties of the dispersion relation 

A plot of k 2 versus a is given in Fig. [| To be noted from the 
figure is that there are always two real values of a, occurring 
in a > and a < — V, respectively. In the region — V < a < 
the dispersion relation returns two real positive values of a, 
provided k 2 exceeds a minimum value, fc^ ln , given by 



4i„ = ^(l + * 1/3 ) 3 



occurring on 

a mi n = -V(l+d 1 / 3 )- 1 > -V 



(15) 



(16) 



For values k 2 < k 2 lin , the solutions are complex conjugate 
solutions in — V < ar < 0. Note that on k = 0, the complex 
solutions are 



«o 



= - V(l + 6)- 1 (1 ± iVS) = a max (l±tVS) 



O 



mm 

-2 




Fig. 2. The dispersion relation for longitudinal waves. Note that 
there is always a real solution in the range a < —V, and another 
real solution in a > 0. In the remaining range — V < a < 
the dispersion relation returns two real values of a only if k 2 



exceeds k^ in . 



For more general values of k 2 < k 2 lin , the general solution to 
the dispersion relation is given through 



aV 1 = — sin 2 4> + iy\ 



with 

k 2 



tair 



k min V tan^ 



and 



tan</.(l + <5 1 / 3 ) 3 



tan </)(v^tan 3 0—1) 
(1 + tan 2 0) 2 (tan0 - v^) 

The range of <j> is restricted to 

5" 1 / 6 <tan0<5- 1 / 2 



(18) 



(19) 



(20) 



(21) 



On tan i 



(5" 1/6 note that k 2 



2 lin and a R = a min , and 
on tan0 = J -1 / 2 , k — and or = a max . Thus the range ( plf ) 
covers the complete wavenumber spectrum where instability 
can occur. In Fig. || we show an enlargement of the dispersion 
relation in the phase velocity range where instability occurs. It 
is obvious that k 2 is basically independent of the phase veloc- 
ity a R , unless a R is very close to a max . Therefore the Taylor 
expansion to first order (Eq. |l0|) is a poor approximation, for 
d J g° R '°^ ~ 0. A second order expansion would provide us only 
with Op which we have already derived in Eq. However, 
Eq. [Kj is very useful for determining the sign of the growth 
rate. Noting that a R < 0, we see that I(a R , 0) is negative for 
(17) negative |& because J|jl > 0. J(a R ,0) is given in Eq. [i~4|, 
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Fig. 3. An enlargement of the dispersion relation for longitudi- 
nal waves in the phase velocity range where instability occurs. 
For small perturbation parameters, 8, k/k m - m is essentially fiat, 
unless or is very close to a max . 



provided a is replaced by hr. Obviously 9J g°^'°' > 
thus a\ < 0. Since k < we get for the growth rate 

7 = kai = \k\ \ai\ > 



< 0, and 
(22) 



As in the non-relativistic case parallel electrostatic waves grow 
as long as the distribution function is inverted, whereas rela- 
tivistic Landau damping would occur forf& >0with P || <0. 
The square of the growth rate is 



7 2 = k 2 V 2 yl = Q- 



(VSt&n 3 (j) - 1)(1 - tsmcfrVS) 



(23) 



(tan 2 </>- l) 2 

In Figs. ^| and ^ we show the growth rate as a function of 
the wavenumber k and phase velocity, or, respectively. The 
growth rate is sharply peaked at wavenumbers close to fc m i n , 
but is essentially independent of the phase velocity, or. 

2.3. Quasi-linear behavior 

The time-dependent behavior of the intensities I(k, t) of the 
excited waves is given by (Lerche [19671 Lee & Ip |l987| ) 

d_ 

dt 



I(k,t) = 2jl(k,t) 



(24) 



Note that instability occurs only for negative wave phase 
speeds or, implying negative wavenumbers fc. 

To describe the long-term influence of the excited waves 
on the beam particles, which are the particles that can resonate 
with the waves, one uses the quasi-linear Fokker-Planck equa- 
tion. Because the longitudinal waves act with an electric vector 
only, and because that vector parallels the magnetic field, the 
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Fig. 4. The growth rate 7 as a function of the wave number k 
for three values of the perturbation parameter 8 = rii/ {rib T 3 ) . 
Most of the instability occurs in a small wavenumber range 
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Fig. 5. The growth rate 7 as a function of the phase velocity or 
for three values of the perturbation parameter 8. The perturba- 
tion parameter controls only the speed of the instability, not its 
spectral form. 



phase space density for the resonant particles then has only its 
momentum parallel to the ambient field influenced by the Ion- 
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gitudinal turbulence. Hence, the corresponding Fokker-Planck 
equation reads 



= JL ( d 2L 

dt dp\\ \ dp\\ 



where the diffusion coefficient, D, is given by 



(25) 



D 



(Apf 

At 



167r 2 e 2 



o 

16ttV I dkI(k)S[k(v { \ -o R )] 

'fcmin 







dk ~\kJ s ^ 



or) 



(26) 



The corresponding diffusion equation for each resonant particle 
species is 

d f Ifi 2 2 9 ( f° „I(k) df\ 

where or is related to the wavenumber k through equation (|l9[). 

Note that some care has to be exercised in treating the 8- 
function in equation (^vj) for two reasons. First, slight varia- 
tions in or (tiii) correspond to large variations in pn for rel- 
ativistic particles. Second, during the evolution of the system 
the initial beam may disperse, for example through interactions 
with electromagnetic turbulence, so that a range of p\\ would 
correspond to a given vu . We consider here two extreme cases: 

(1) when the variation in phase speed is much larger than 
the variation in transverse speed; 

(2) when the variation in transverse speed is much larger 
than the variation in phase speed. 

Consider each case in turn. 



2.4. Highly variable phase speed 

In this situation one can effectively set the resonant particles to 
a cold gas approximation and deal only with the bulk flow. 

One conventionally argues that the wave generation is suffi- 
ciently rapid that the local (in time) wave intensity has reached 
a quasi-steady state in the sense that it alters quasi-adiabatically 
as the particles change their bulk streaming parameters to feed 
the waves. 

Because of the unknown nature of the initial wave spec- 
trum, and because other processes ignored in the development 
will also influence the generation of waves, one of the standard 
devices is to ignore the wave intensity spectrum generation and 
to use models of how one believes the waves spectrum has 
evolved to its current state. This particular decoupling of the 
wave and particle behaviours has been an enormous success in 
cosmic ray astrophysics in general and in heliospheric physics 
(A thorough review can be found in Schlickeiser, 2001 ). Under 
this scenario one writes equation ( j26[ ) now with I(k, t) a wave 
intensity spectrum to be prescribed by the user. 

In addition, the parameters related to the relaxation of the 
plasma, V, or, etc., are considered to be independent of time in 
evaluating the diffusion coefficient (26), but are then allowed 
to evolve with time once the diffusion coefficient is evaluated 
in order that the long term evolution of the plasma distribution 
function can be tracked. 



The particle velocities, vu, for which resonance can be es- 
tablished, are those which satisfy 



-V 



< v\\ < 



-V 



1 + 5 ~ " 1 + S 1 / 3 
corresponding to particle momenta, pu, in the range 



-VT 



< P -l< 



-VT 



(28) 



(29) 



V /1 + <5T 2 (2 + <J) " m ~ + 5 1 ' 3 T 2 {2 + 8 1 / 3 ) 

Clearly, the range of resonant particle momenta is sensitive to 
2ST 2 and 25 1 ^ 3 T 2 . If both are small compared to unity then the 
range of resonant momenta is close to the initial momentum of 
the beam particles, — P, and has a width of ~ P5 X / 3 T 2 , while, 
if 25T 2 and 25 1 ^ 3 T 2 are both large compared to unity then the 
particle momentum range is far from — P and has a width of 
~ Pj V26T 2 . Thus the quasi-linear evolution is sensitive to 
the values of 25T 2 and 25^ 3 T 2 , both of which are required to 
be small compared to unity in the astrophysical application to 
AGN (Pohl & Schlickeiser |2000| ). 

When we consider the growth rate (Eq.^3[) as a function 
of phase velocity, we see in Fig. || that the growth rate is es- 
sentially independent of or. Neglecting the possible effects of 
damping and cascading, we can therefore expect that the inten- 
sity spectrum in phase velocity, /(or) = I{k)-g£- is flat be- 
tween a m j n and a max and zero outside this range (Akhiezer et 



al. 1975). Upon interaction with this wave spectrum the beam 



will relax to a plateau distribution between the velocities — V 
and vu = a m i n . The energy available for the build-up of the 
waves comes from the beam. Because the wave growth is much 
faster than the relaxation of the beam, we may estimate the final 
wave intensity spectrum as the ratio of the energy lost by the 
beam during relaxation to a plateau between the velocities — V 
and v\\ — a m in and the phase velocity range in which effective 
growth occurs, namely S 1 ^ 3 V. Therefore we may write 

2 p 
Aa R ~ V6 1/3 

AE n ? m„c 2 r 3 

/(OR) ~ ~ T~TZ = / 



Aa R 



2V 







(30) 
(31) 
(32) 



where we keep in mind that ajx is to be taken between a m j n and 
B max . The diffusion coefficient is then 



D(p«) = 16ttV/, 



daR W\ 6 ~ ° R ' ) 



= 167r 2 e 2 



IMflR = v \\)\ 



2tt m e m p c 2 T 3 — 
n b 



(33) 



where we have used again that k ~ k m - m ~ — V/Q in the 
resonant range. The time scale for the relaxation of the beam 
by the electrostatic instability can then be estimated as 



p2 

D 



m p rib 



2n Top m Q r 



5 • 10" 



Yrii 



sec 



(34) 
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In term of the interstellar matter density in the laboratory frame, 
n*, and the scaled variables rig = 10 _8 nh and T2 = 10~ 2 T the 
relaxation time scale reads 



5- 10" 



T 2 n* 



sec 



(35) 



which is two orders of magnitude shorter than the time scale for 
the electromagnetic instability derived in Pohl & Schlickeiser 
(120001). 



2.5. Highly variable transverse speed 

Suppose the initial beam is slightly dispersed, such that all par- 
ticles have the same total momentum, P, but instead of flying 
exactly backwards (/i = — 1) the particles homogeneously oc- 
cupy a fixed solid angle element defined by [i < fi c . This initial 
condition approximates the situation after the onset of pitch an- 
gle scattering but long before the isotropization is completed. 
The initial distribution function in the blast wave frame is then 



f(p,t = 0) = 



rihS(p±)5(p 



2iTp± 



pf,-p)eu c 



2ttP 2 (1- 

The dispersion relation (||) becomes 



n 2 



ST 2 



dy 



<*(1 + Mc) ./-l (y-%) 



(36) 



(37) 



where (3 — V/c,y = P\\ / P, and 5 — m / (rib T 3 ) as in Eq.(|l7j). 
The integral in Eq.([37|) can be trivially solved and we obtain 
the dispersion relation for a dispersed beam with x — a/V and 

e = 1 + fj, c 



k 2 V 2 

n 2 

with 

SO) : 



1 



8T 2 g{x) 



(38) 



dy 



(y-vY 



which for e — > would reduce to Eq. |14| For small e we can 
derive modifications 



x 



x + Xl 



(39) 



to the cold beam solutions x (Eqjl8|) of the dispersion relation 
by expanding the right-hand-side of Eq.B8l RHS(a;, e). 



k 2 V 2 aRHS(x ,e) 
— e- 



Q 2 



e=0 



de 

x 2 a 2 RHS(x,0) 



X\- 



<9RHS(:r,0) 



dx 2 



dx 



RHS(a;o,0) 



(40) 



The second order expansion in x is necessary because the first 
order term in x is zero for tan 6 = <5 -1 / 6 . 
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Fig. 6. The normalized phase velocity modification, Sftxi, for 
different beam openings, e, and different bulk Lorentz factors, 
r. The perturbation parameter 5 is held constant. The phase ve- 
locity modification scales approximately linearly with e, unless 
aR is close to 

The term RHS(xo,0) is identical to the right-hand-side 
of the dispersion relation for an ideal beam (Jl4t), and hence 
equal to k 2 V 2 /£l 2 . The remaining binomial expression in x\ is 
solved by 



xi = (1 + x ) 



(1 + xy 1 ) 3 + 6 - ^/fj^ej 



3[(1 



-1^4 



(41) 



where 

f(x ,e)= [(l + x^f + S] 2 

1 + x 
Here we have also used 



lim xi =0 

to resolve the ambiguity in sign. 
For the expansion we require 

< 1 



1+x 



(42) 



(43) 



(44) 



which leads to a limit for the opening half-angle of the beam 

9 = — 7r + arccos fj, c 



e ~ 9 2 < <5 1/4 tan-50 = 5 1/2 for tan0 = ^ 1/2 



(45) 



The limiting beam opening half-angle is about 10 arcminutes 
for S = 10~ 10 . For all e satisfying Eq||5]the parallel velocities 
v» are always larger than the phase velocities in the range where 
instability occurs. 
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Growth rate modification 
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Fig. 7. The growth rate modification, 71, as a function of the un- 
modified phase velocity ar for different beam openings, e, and 
different bulk Lorentz factors, T. The perturbation parameter 
S is held constant. The growth rate modification has a similar 
dependence on the parameters as has the phase velocity modi- 
fication. 
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Fig. 8. The modified growth rate, 7 + 71, as a function of the 
modified phase velocity, <2r + V$txi, for different beam open- 
ings, e. The perturbation parameter S and the bulk Lorentz fac- 
tor r = 300 are held constant. Note that this plot has a linear 
scaling in units of 10 -4 . The solid line applies for all e < 10~ 8 , 
for which the modification has no visible effect. The dotted line 
represents the enhanced growth rate for e = 10 -6 . 



In Fig.g we show the phase velocity modification, V^Rxi, 
for different parameters. The phase velocity modification scales 
approximately linearly with e, unless or is close to a m j n , which 
reflects the dominance of the first order term in x\ over the sec- 
ond order term in Eq.^ for all not too close to a mm . Fig^ 
shows the growth rate modification, 71, for different parame- 
ters. 

In Fig.||we show the modified growth rate, 7 + 71, as a 
function of the modified phase velocity, <2r + VSlxi. To be 
noted from the figure are 

- The growth rate increases if the beam has a finite opening 
angle. 

- The phase velocity range where instability occurs is only 
marginally extended to phase velocities smaller than a m i n . 

3. Discussion 

We have calculated the excitation of longitudinal, electro- 
static turbulence subsequent upon the injection of a relativistic 
electron-proton beam into a cold background plasma. In con- 
trast to electromagnetic waves, which scatter the particles in 
pitch angle but preserve their kinetic energy until the distribu- 
tion is isotropized, the electrostatic waves change the particles' 
energy until a plateau distribution is established. Our results 
can be applied to models of AGN and GRB, which are based 
on a relativistic blast waves traversing the ambient medium. 

As one example of such models, two of us have recently 
published a calculation of an electromagnetic two-stream insta- 
bility in a channeled relativistic outflow, a jet, and shown that 
the isotropization of the incoming interstellar particles would 
provide relativistic particles in the AGN jets, whose high en- 
ergy emission has characte ristics typical of BL Lacertae ob- 
jects (Pohl & Schlickeiser 200C ). The calculations presented 
here expand on these results. We have shown that for param- 
eters suitable to explain the high energy emission spectra of 
AGN, the electrostatic instability is much faster than the elec- 
tromagnetic instability, even in case of a weak dispersion of 
the incoming beam. The beam will therefore relax to a plateau- 
distribution in p||, which is, however, still unstable with re- 
spect to electromagnetic waves. The spectrum of electromag- 
netic waves excited by a plateau-distribution in p\\ will not be 
a ~kr 2 spectrum as in the case of a cold beam, for it can be 
understood as a superposition of the wave spectra produced by 
cold beams with different energies, each of which would be a 
k~ 2 spectrum with energy-dependent upper and lower limits 
in k Nevertheless, isotropization in the collimated blast wave 
would still occur on roughly the same time scale as calculated 
by Pohl & Schlickeiser ( |2000| ). 

Therefore, the main effect of the electrostatic instability is 
a change in the spectrum of swept-up particles. Instead of the 
prior result for the differential sweep-up rate in the absence of 
the electrostatic instability 
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Fig. 9. 7T° -decay spectra calculated with and without the elec- 
trostatic instability. The spectra refer to a observed time of 
0.125 hours after the blast wave ran over an isolated cloud. At 
this time the spectrum of swept-up particles has evolved very 
little. At the peak, about 50% of the beam particle energy is 
lost due to the effect of the electrostatic instability. There is 
little effect at smaller photon energies. 



In the radiation spectra this change corresponds to a reduc- 
tion of efficiency by roughly 50% at the highest energies and 
little change at smaller energies, because in the inelastic pp 
collisions the pion source power scales roughly with the ki- 
netic energy of the incident proton, hence the pion spectra are 
dominated by the pions produced by the protons of high en- 
ergy for flat proton spectra. In Fig.^ we show the spectra of 
7T° -decay 7-rays for the case of a fast electrostatic instability 
(Eq.^7|) in comparison with the case of the electromagnetic in- 
stability only (Eq.|46|). At photon energies well below the peak 
energy of the vF v -spectra the difference in spectral index is 
only As ~ 0.1. Similarly, small changes can be expected for 
the injection spectra of secondary electrons and hence for the 
synchrotron spectra at optical to X-ray frequencies. Therefore 
the isotropization and radiation conclusions drawn by Pohl & 
Schlickeiser (200C ) are still valid. About half of the energy flux 
carried by the incoming beam is channeled into electrostatic 
turbulence. It is instructive to calculate the group velocity of 
the waves with which the wave energy is transported. In Fig . [U)| 
we show the group velocity for the initial ideal beam (Eqjl9|) 
in the phase velocity range where instability occurs. The group 
velocity is fairly independent of the weighted beam intensity, 6, 
and it is always larger than 0.5 V ~ 0.5 c. Therefore, the energy 
carried by the electrostatic waves is rapidly transported toward 
the backside of the blast wave. Because the waves carry 50% of 
the energy flux of the incoming beam with a transport velocity 
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Fig. 10. The group velocity of electrostatic waves as a function 
of phase velocity for the initial ideal beam. The group velocity 
is always between about 50% and 100% of the beam velocity, 
which result is fairly independent of the weighted beam inten- 
sity, 5. 



of 50%-100% of the beam velocity, the energy densities in the 
final wave field and the initial beam should be similar. 



Uw = J dan I(clr) — Hi m p c 2 T = n* m p c 2 T 2 



(48) 



The high energy density of waves has two important aspects, 
which are briefly discussed in the following sections. 

3. 1. Nonlinear evolution of the wave field 

Quasi-linear theory is not valid under all physical conditions. 
When the energy density of the wave field is not small, nonlin- 
ear effects may become important. 

One important effect is nonlinear Landau damping 



(Sagdeev & Galeev 1969 ). This is the interaction between ther- 
mal particles and the beat of two electrostatic waves: if the 
phase velocity of the beat wave is almost equal to the thermal 
velocity of charged particles in the background plasma, the par- 
ticles would Landau-damp the beat wave and hence be heated. 
If the remaining wave has a superluminal phase velocity, it 
will not resonate with the beam particles anymore. As a conse- 
quence, the number of scattering processes between the beam 
particles and the waves, which are responsible for the produc- 
tion of a plateau in the beam distribution function, would be 
drastically reduced. This process reduces the wave intensity in 
the resonant phase velocity range as long as the intensity of 
waves with superluminal phase velocity is not so strong that 
nonlinear interactions in this regime cause wave energy to be 
retransferred to resonant subluminal phase velocities. 

The nonlinear Landau damping is a process resulting from 
weak turbulence theory: although waves interact, the wave en- 
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ergy is small enough that individual waves are identifiable. If 
the wave energy is too high, strong turbulence theory applies 
(Zakharov 1972). The charge distribution in the background 
plasma, which results from the electric field of the waves (es- 
sentially the RHS of Eq.Q), becomes so strong that it affects 
the dispersion relation. As a result, the wave field will collapse 
into cavitons (Goldman 1984) in which the plasma has a low 
density. This effect will depend on the temperature of the back- 
ground plasma, and on whether or not the diffusive particle mo- 
tion can effectively smear out the depleted cavitons, and thus 
impede their build-up. 

Without having performed a detailed study of these effects 
in a background medium with possibly time-variable tempera- 
ture caused by heating processes, we can only speculate on the 
impact of the nonlinear effects, which would increase the time 
scale for the plateau production of the beam (Eq.[35]). However, 
because that time scale is more than two orders of magnitude 
shorter than the time scale for the electromagnetic instability, 
even if important such nonlinear effects would likely still per- 
mit the production of a plateau on a timescale much faster than 
the isotropization. 



The existence of radio galaxies implies that the collimation 
of AGN jets is well maintained on linear scales of hundreds of 
kiloparsecs. A semi-relativistic plasma jet could hardly remain 
collimated on these scales, because the expansion by internal 
pressure alone would cause a substantial opening of the jet. We 
may take the observed collimation of jets in radio galaxies as 
empirical evidence against excessive heating. For the fate of 
the electrostatic wave field, this observation indicates that the 
reflection efficiency of waves at the backside of the blast wave 
is small, whatever the cause. 

Since the waves carry momentum, even a small reflection 
would transfer momentum to the background plasma at the 
backside of the blast wave, which would then "boil off the 
blast wave and fill its wake. The waves that have escaped the 
blast wave would fill the wake and the regions around it, and 
be finally damped by the interstellar plasma in the host galaxy 
of the AGN, thus giving rise to additional heating of the inter- 
stellar medium around the AGN jet. 
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3.2. Where does the wave energy go? 

Part of the wave energy will be transferred to the background 
plasma in the blast wave. We have already discussed the nonlin- 
ear Landau damping, in which thermal particles would absorb 
the beat of two electrostatic waves. But, in the general case, the 
natural width of the resonance Lorentz profile combined with 
the Maxwellian distribution function allows thermal particles 
to Landau damp the waves with v$ ~ c, albeit with a small 
probability. 

We do not know the fraction of the wave energy density that 
is channeled into heating processes of the blast wave plasma. 
Here we concentrate on a qualitative discussion of the possible 
effects. 

The wave energy is transported quickly toward the backside 
of the collimated blast wave. At the transition between the blast 
wave and its wake, the dispersive properties of the blast wave 
change. We can, therefore, expect that a part of the wave field 
is reflected at the backside, and possibly again reflected at the 
front side and so forth. If the reflection coefficient is not small, 
the wave energy density inside the blast wave would increased, 
until nonlinear effects at the front and backside cause the waves 
to escape. The heating rate would be similarly enhanced. 

What would be the effect of strong heating in the blast 
wave? As an extreme example consider that the full wave 
energy flux is available for homogeneously heating the blast 
wave. The heating rate then is 

^10 6 ^j K/sec (49) 

where we have used the scaled variables ng = lCP 8 rt&, d\j, = 
10~ 13 d, and T 2 = 10~ 2 r. This heating rate would be too 
large to be balanced by any cooling process of a non-relativistic 
plasma. Therefore run-away heating would occur until non- 
thermal cooling processes set in at semi-relativistic energies. 
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